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Abstract 

We investigate the evolution of the Universe filled with barotropic perfect fluid in Eddington- 
inspired Born-Infeld gravity. We consider both the isotropic and the anisotropic universe. 
At the early stage when the energy density is high, the evolution is modified considerably 
compared with that in general relativity. For the equation-of-state parameter w > 0, the 
initial singularity is not accompanied as it was discovered for radiation in earlier work. More 
interestingly, for pressureless dust {w = 0), the initial state approaches a de Sitter state. 
This fact opens a new possibility of singularity-free nature of the theory. The anisotropy is 
mild, and does not develop curvature singularities in spacetime contrary to general relativity. 
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1 Introduction 



Einstein's theory of general relativity (GR) formulated in 1916 is very successful in agreement 
with many phenomenological and experimental results. However, it is well-known that GR suffers 
from the singularity problem which seems to be unavoidable in the beginning of Big Bang, or at 
the center of black holes. Very recently, Banados and Ferreira suggested an alternative theory 
inspired by Eddington's theory of gravity [1]. This Eddington- inspired Born-Infeld (EiBI) theory 
of gravity requires only one more parameter k other than the gravitational constant G, which is 
reviewed below. 

In Ref. [T], the authors showed that the EiBI theory in vacuum is equivalent to GR, while it 
deviates from GR in the presence of matter. Most interestingly, the Universe driven by radiation 
is free from the initial singularity; the Universe experiences a bouncing with a finite size for k < 0, 
or there is a state of minimum size for which one takes infinite time to reach from the present for 
K > 0. The latter is interpreted as the "nonsingular initial state" of the Universe. 

In Ref. [2], the authors considered the modification of Poisson equation in EiBI gravity, and 
obtained singularity-free solutions for the compact stars composed of pressureless dust and poly- 
tropic fluids. In Refs. [Sill], the cosmological and astrophysical constraints on the EiBI theory 
was studied. In Ref. [5], the constraint on the value of the coupling parameter /t was investigated 
by using the solar model; the result does not rule out the EiBI theory as a possible alternative to 
GR. A number of subsequent articles studied the tensor perturbation [6] , bouncing cosmology [7] , 
the five dimensional brane model [S], the effective stress tensor and energy conditions [H] in EiBI 
theory, etc. 

The EiBI action considered in Ref. [1] is given by 

S'EiBI = ^y"^'^[ \l-\9i^v + ^^M-(r)| - \\|-\g^.^ ] (1.1) 

where \g^v\ denotes the determinant of g^y, A is a dimensionless parameter which is related with 
the cosmological constant, and SvrG was set to unity. Then this theory becomes a one-parameter 
(k) theory. In this theory the metric gfj_u and the connection F^^ are treated as independent fields 
(Palatini formalism)]^ The Ricci tensor R^^iV) is evaluated solely by the connection, and the 
matter filed $ is coupled only to the gravitational field g^^^. 

According to the Palatini formalism, one should consider the equations of motion by varying 
the action (11. ip with respect to (w.r.t) the fields g^y and F^^^ individually. Variation of the action 



^In the original Palatini formalism |10| . the matter action Sm depends on F^^, as well as g^i,^ and the connection 
is not symmetric (there exists a torsion Tj'^j^j). However, in EiBI theory, Sm is assumed to depend only on 5^1/ and 
the torsion is assumed to be absent. 
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w.r.t. gfj_i^ leads to the equation of motion, 

-[{g + kR)-']^"' - \g>"' = -kT>"', (1.2) 



where [{g + nR)~^Y'^ denotes the matrix inverse. The energy-momentum tensor T^'^ is given by 
the usual sense, 

T'^^ = (1.3) 
^J-\g\ og^u 

For the variation of the action w.r.t. F, one introduces an auxiliary metric g^^ defined by 

Qfiu = g^u + K.R^lu■ (1-4) 
Then the variation of the action (11.11) w.r.t. the connection F^^ gives 

Kr = 0, (1.5) 

where g'"^ = (g~^)'"^ is the matrix inverse of qpa-, and V'" denotes the covariant derivative defined 
by the connection F. This equation is the metric compatibility which yields 

r^/J = ]^q^"{qaa,p + qi5a,a + Qa^^a)- (1-6) 

Therefore, Eq. (11. 4p can be regarded as the equation of motion since the Ricci tensor is evaluated 
in terms of qpa through the relation (II. 6p . Using Eq. (II. 4p . the first equation of motion (II. 2p can 
also be simplified. 

We would like to mention a couple of properties of the equation (II. 7p . First, when T'^'^ = 0, the 
metric satisfies the relation g^^ = g^^/A. Then Eq. (II. 4p becomes i?^,^ = Ag^,,, where A = (A— 1)/k. 
This implies that the EiBI theory reduces simply to GR in vacuum. Second, the matter field in 
EiBI couples only with the metric (7^^, so the conservation law V^T'"^ = is expected to hold. 
Here, denotes the covariant derivative defined by the Christoffel symbol based on g^^iy. In 
Appendix A, we show that this really holds from Eq. (II. 7p . 

In this paper, we investigate the Universe filled with perfect fluid in EiBI theory. The per- 
fect fluid drives the Universe in a different manner from that in GR, since the effective energy- 
momentum tensor is differentia We precisely investigate the evolution of the Universe case by case 
depending on the equation-of-state parameter w = P/p for barotropic fluid. We also investigate 
the Kasner-type anisotropic universe. We analyze differences from as well as similarities to GR in 
the results. 



q'^'^ = A^^^ - kT"". (1.7) 



^In Ref. [9], the effective energy- momentum tensor of perfect fluid in EiBI was studied. However, the geometry 
part responding to the effective energy-momentum tensor was described by the auxiliary metric. Therefore, the 
evolution of the Universe was not very evident. 
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2 Field Equations with Perfect Fluid 



In this work, we consider barotropic perfect fluid of wliicli tlie energy-momentum tensor is given 
by 

T^'^ = (p + + pg''\ (2.8) 

For the Kasner-type anisotropic universe, the general metric ansatz can be 



and the auxiliary metric can be 



Qf^^dx^dx" = -X^dt^ + Y 



(2.9) 



(2.10) 



where fl, /3±, /3±, X, and Y are functions of t only. 

With the above metrics, the nonvanishing components of the equation of motion (11.71) are 



XY 



Y^ 

e^^X 
X 



+ A 



g3n+2(/3++v/3^_) g2n+2(/3++V3/3-) 

XY A 



g3n+2(/3+--V3/3_) g2t7+2(/3+-V3/3_) 



g2Q+2(^++x/3/3_) ' 

up 

g2C+2(/3+-V3/3_) ' 



g3r2-4/3+ g2r2-4/3+ 

From Eqs. fl2A2|) and fl27[3|) . we get 



g2n-4/3+ ' 



/3_ = /3_ 



and 



(A-«:p)e2(^+-/5+). 



Plugging these relations into Eq. (I2.14p . we get 

/3+ = and XY = {X - Kp)e^. 

From Eqs. (12. lip and (I2.16p . we have 

X = i4^1^^, and Y =[{X-Kp){X + Kp)f^e''. 



(2.11) 
(2.12) 
(2.13) 
(2.14) 

(2.15) 
(2.16) 



(2.17) 



(A + «:p)V4' 

With Eqs. ( I2.15P and (12.160 . the nonvanishing components of the equation of motion (11.40 
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become 



- + 1 = 3k 



y2 _ g2n 



y2 _ g2n 



d [y 

Jt \ Y 

dt [y 

dt I Y 



Y 
Y 



XY 
XY 



2(/3U/3- 



3— - — 



3^ 


X 




^ X 




X 




^ X 




-2/3. 





(2.18) 

(2.19) 
.(2.20) 
(2.21) 



Equating fl2.19l) - fl2.2ip for the anisotropic factor, we get 

: 



X 



(2.22) 



where c± is an integration constant. For the isotropic expansion, c± = 0. Manipulating Eqs. f l2.18p - 
fl2.2ip with Eqs. f l2.17p and (12. 22 p . we get two equations of motions for X and Y, 





1 




6k, 



d Y 

where = ci + and 



1 + 2X' 



XY 1 
XY ~ 2^ 



3X^ 



(A-Kp)2j ye 



/(x,y,p) 



6k 



X' 



{x-Kpy 



2^' 



(2.23) 
(2.24) 



/(X,r,p) = 1 + 2X2 



3X^ 2X2 

(A - Kpy + ^'"'^ F^- 



(2.25) 



As it was mentioned earher, the conservation law for the matter is given by V^T^'' = (see 
Appendix A). For perfect fluid, it reduces to 



p + 3n{p + p) = 0. 



(2.26) 



For barotropic fluid, the equation of state is given by p = wp, and the solution to the above 
equation is given by 

p = poe-=^(i+-)^ = poa"'^'+"'\ (2.27) 
where we defined a = which we shall call the scale factor. 
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Using Eqs. (KT7} . ([226]), (KTf} . and p = wp, the Friedmann equation of first kind (K23\i 
becomes 

^ ^2^(A^^ (2.28) 



X 



6k 

(A + Kpf + 2(A - WKpf/\X + Kp)'^/^ - 3(A - WKp){X + np) + QKc^{p/pof/^^+'"^ 
[{3/A)kw{1 + w){X + Kp)p - (3/4)fi;(l + w){X- WKp)p + (A - WKp){X + Kp)f 



and the equation for the anisotropic part fl2.22p becomes 

3 Evolution of Universe 

In this section, we investigate the evolution of the Universe for various values of the equation-of- 
state parameter w by analyzing the two field equations f l2.28p and (12.291) . We shall focus on the 
case of K > 0. The parameter A is related with the cosmological constant. Although the value 
of A is not restricted, in order to see the pure role of perfect fiuid, one can set the cosmological 
constant to zero (A = 1). We discuss our results mainly for A > 0. 

3.1 w>0 

3.1.1 Nonsingular Initial State 

When X — wnp = 0, i.e., at p = ps = X/wn, the expansion rate in Eq. (I2.28P becomes zero, H = 0. 
Let us expand about this point. We can write the energy density and the scale factor as 

p = Pb — £, and a = as + e, (3.30) 

where e and e are small quantities. From Eq. (I2.27p . we have 



(3.31) 



From this, we get the relations, 

pB = poO'B^^^^\ and e = 3(1 + w)pBa~^^e. (3.32) 
With the aid of these relations, in Eq. (12.280 can be expanded as 



, 2 

a 



, 8kW^ Ul+wfX^ + 6KW^c\pB Pof/^^+''^ 0,^, . . 

^ ~ .... =Hl{ 1 , 3.33 

27 1 + ti; 4A4 ^\aB ' ^ ' 
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where „ „ 

^0 = 7^ + TT^- 3.34 

Note that the first term in Hq comes from the EiBI correction, and that the second term comes 
from the anisotropy with EiBI correction. From Eq. (I3.33p . we finally obtain the scale factor 

a^aB + Ae^-^«*, (3.35) 

where A is an integration constant. We note that the expanding solution is possible only for A > 
considering the definition of X and Y in Eq. f l2.17p . Therefore, we have the expanding solution 

a{t) ^aB + Ae^"*. (3.36) 

As it was studied for radiation {w = 1/3) in Ref [T], this solution indicates that there could exist 
a nonsingular initial state of the Universe. The Universe begins with a finite scale factor a^, for 
which the Universe has a maximum value of energy density pb = X/wn. However, it takes infinite 
cosmological time to reach this state flushing back in time. Therefore, there would be no horizon 
problem. From our result, this is true not only for the radiation-dominated universe, but also for 
all the cases of > if k > 0. (See Fig. 1.) 

The Universe undergoes accelerating expansion in the beginning, even when the cosmological 
constant is absent (A = 1). Later at the low-energy regime, the Universe expands with deceleration 
as we shall see in the next subsection. The e-folding depends on the parameters involved, but 
from numerical calculations we observe that it is order of unity ~ 0{1) as a whole. The e-folding 
becomes considerably large ~ (9(10) as — 0. This can be interpreted as a limit of dust [w = 0) 
in the next subsection. 



3.1.2 Late-Time Evolution (p < X/wk) 

When the energy density of the Universe becomes small as the Universe expands, the Friedmann 
equation f l2.28p approximates as 

The first term corresponds to the cosmological constant, the second term is the linear dependence 
in p similar to GR (but for X ^ 1 there is an EiBI correction in the coefficient), and the third 
term is the correction purely from the anisotropic expansion. The anisotropic term is dominant 
over the linear term for w > 1 in the low-energy limit. 



X-1 1 {X-l)w{w + l] 
3k ^3 2A 
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When the cosmological constant is absent (A = 1), the expansion becomes 



2 

1 / n \ 



^^-3"+^^; . (3,38) 

which has no K-dependence. Therefore, we can conclude that the late-time expansion of the 
Universe approximates to that in GR. When the cosmological constant is present (A 7^ 1), it must 
dominate the late-time low-energy universe. Therefore in the end, for A > 1 the Universe must 
asymptotes to de Sitter, and for A < 1 to anti-de Sitter. 

Let us briefly discuss the case of A < 1. As one can see from Eq. f l2.28p . or fl3.37p . there exists 
a critical density p = Pb ^ Pb = \/wk at which the expansion stops, H'^{pb) = 0. After the 
moment of p = pb, < and the Universe contracts. Near pf,, the Hubble parameter behaves as 

(X p — pb (X ttb — a, (3.39) 

and the scale factor becomes 

a = ab- HbiU - tf. (3.40) 

Therefore, aX t = % the expanding universe bounces back to contract. This type of late-time 
behavior is more or less similar for other values of w, except w < —1 ioi which the late-time 
universe corresponds to the high-energy state. 

3.2 w = (Dust) 

When p = 0, the expansion behavior is very peculiar. We precisely analyze this dust-filled universe. 
In this case, the expansion rate fl2.28p becomes simpler, 

, 8 -2A^ - K\p + 2A^/^(A + Kpfl^ + (^^ + Qnc^/pDp' 
^ = 3 .(.p + 4A)^ • ^^-^^^ 

There is no = point for A > 1, which is different from the > case. Note that the 
numerator can be reexpressed as 



A2^2 

where S = y^l + /tp/A > 1. The right-hand side is positive definite if A > 1, so > 0. (For the 
isotropic case, the scale factor a can be obtained explicitly. Please see Appendix B.) 



(^-1)^ + 2(A + 1)(^-1) + 2(A-1) +6fic'{-^j >2\^S\\-1) + 6kc' i^-^] ,(3.42) 
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3.2.1 Early-Time Evolution (p > X/n) 

In the high-energy limit, the expansion rate in Eq. (13.41 p becomes 

1 



1 & 
6 \ K k'^Pq 



16 /Ay/' 



.1/2 



(3.43) 



Very interestingly, in the high-energy limit, the Universe approaches the de Sitter state (see 
Fig. 2.), 

Acff 

3 ' — ' ~ ' ^ ~ 



(3.44) 



The effective cosmological constant originates from the EiBl nature of dust with contributions 
coming from isotropy as well as anisotropy. This means that for dust in high density, the repulsive 
gravity is produced in EiBI theory. This provides a new interesting scope of singularity-free nature 
in EiBI theory. First, the initial singularity is not accompanied in the dust-filled universe. This 
singularity-free initial state is somewhat different from what was obtained for the w > case; it 
is de Sitter state. Second, this repulsive nature of gravity suggests a new possibility of avoiding 
the singularity formation in collapsing dust. When pressureless dust collapses gravitationally and 
reaches the high-density regime, the repulsive nature of EiBI gravity may arise to prevent further 
collapse. 

3.2.2 Late-Time Evolution (p < X/k) 



In the low-energy limit, the expansion rate becomes 



X-1 
3k 



k(3 - A) 
16A2 



A2p2 



Poi 



p' + 0{p' 



(3.45) 



Therefore, the late-time expansion is similar to that in GR. 



3.3 



1/3 < w; < 



For this case, there exists a moment at which H becomes singular. The denominator of in 
Eq. fl2:28|) vanishes at 



P = Pc 



A (1 - w)(l - 3w) ± (1 + w)Vl-4:2w + 9w^ 



(3.46) 



-Aw{l + 3w) 

For the negative root, there is no singular point since pc is negative for all values of w in the range 
if A/k > 0. However, for the positive root, pc is positive definite. At this value of the critical 
density, the H is divergent. (See Fig. 3.) The Hubble parameter around pc takes the form. 



[a — ttr 



(3.47) 
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where Oc is the scale factor at p = pc and h is a constant determined from Eq. f l2.28p . The scale 
factor is solved as 



a{t) ^ac± ^/2a^h\t-t^\, (3.48) 

where tc is the time of the critical moment. The scale factor is finite, a = ac, at the critical moment, 
but the expansion rate H diverges. Therefore, a curvature singularity is formed at that critical 
moment. The Universe is divided into two sectors by this critical moment p = p^. The former 
high-density universe ends up with the singularity within finite time, and the latter low- density 
universe begins with the singularity. 

3.3.1 High-Density Universe (p > pc) 

In the high-energy limit p ^ \\/wn\ of this high-density universe , the expansion rate in Eq. (12.281) 
becomes 



1^ \Po 



2w 



3/2 „, "J*- { P \ 



O^p""). (3.49) 



3(1 + 3u;)2 

For — l/3<w<0, the first term which is linear in p, is dominant and the evolution is similar to 
that in GR. (The coefficient is a bit different.) The expansion is power- law, 

/ ^ \ 3(TT^ 1 + 3w 

a{t) ^ - , where to = -====——-. (3.50) 

V'^o/ V3(-w)^/^Po(l + w) 

The expansion power is 2/3 < 2/3(1 + u;) < 1, so the Universe undergoes decelerating expansion. 
As it was mentioned, the Universe will end up with a singularity at p = pc while it approaches a 
finite size. 

3.3.2 Low-Density Universe (p < pc) 

The low-density universe begins with a singularity at p = pc from a finite size. At late times in 
the low-energy limit p <^ \X/wk,\, the Universe approximates to that in GR, 



3ac 



1 _ {\-l)w{w + l) 
3 2A 



p + 0(/). (3.51) 



In the absence of the cosmological constant (A = 1), the Universe decelerates. 
3.4 -l<w;<-l/3 

For this case, there is no singularity in since pc in Eq. fl3.46p is negative. In the high-energy 
limit, the expansion rate is in the same form fl3.49p . (See Fig. 4.) For this case, however, the 
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a(t) 



second term becomes dominant which is the anisotropic correction from EiBI. The scale factor is 
given by 

t\~^ , (1 + 3w)k 

- , where = ^-^ -■ 3.52 

The expansion power is 1/3 < —l/3w < 1, so the Universe undergoes decelerating expansion. For 
the isotropic case, the second term is absent and the evolution is similar to that in GR. 

At late times in the low-energy limit, the expansion rate is again in the same form fl3.5ip . 
The evolution of the Universe approximates to that in GR. In the absence of the cosmological 
constant (A = 1), the Universe accelerates in this case. 



3.5 w = —1 

For this case, the perfect fluid corresponds to the cosmological constant. The expansion rate is 
given by 

r2 _ ^ + ^ + I «-6^ 

3k 3 (A + upoy 



(3.53) 



which is exactly the same form as in GR. (Although the EiBI parameter k appears in the 
anisotropic contribution, it can be absorbed since the integration constant c± is arbitrary.) 
For the isotropic case, the expansion is exponential as usual, 

A - 1 



a{t) = ev"^/^*, where A 
For the anisotropic case, the scale factor is given by 

, , r 3c2 

a(t) = ^ 

and the anisotropic part is also obtained exactly, 



Po- 







sinh^/^ ( 


\/3At) 



e^±W = tanh" 



'3A 



(3.54) 



(3.55) 



(3.56) 



3.6 w; < -1 

This case corresponds to the phantom matter in GR. The solution (I2.27P to the conservation 
equation tells that the energy density increases as the universe expands, i.e., as a = increases. 



(3.57) 



Therefore, the low-energy limit corresponds to the early universe. 
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3.6.1 Early-Time Evolution (p < 

At early times, the energy density is low. The expansion rate becomes 



1 

A2 



3w )k, 



2A2 



P + 0{p' 



2 

1+w 



A- 1 
3k 



(A — l)w{w 



2A 



(3.58) 

For the isotropic case (c = 0), the expansion is very similar to that in GR. (See Fig. 5.) When there 
is an anisotropic expansion, the first term (c^/ A^-term) is dominant which is inversely proportional 
to the energy density. The second dominant term is the second term for —3 < w < —1, and is 
the cosmological constant term for w < —3. If we consider the most dominant term only, the 
expansion rate becomes 

i/2 ^ ( il ) = l^a-\ (3.59) 



(3.60) 



y \poj y 

which behaves like the stiff matter in GR, and the expansion at early times becomes 



a{t) 



{t - to) 



where to is an integration constant 



i 



3.6.2 Late-Time Evolution (p > \X/wk\) 

At late times, the energy density is high. The expansion rate becomes 



3(1 + 311^)2 



P + 



3(1 + 3w)k 



2{-w)-^/^{3w^ - 10-»7-9)A 
3{l + 3wy^K 



o{p-'). 



(3.61) 



The anisotropic contribution is negligible, and the expansion is similar to that in GR. When only 
the first term which is most dominant is considered, the expansion at late times becomes 



a{t) 



tr-t \ 



to 



r<0 



where 



tn 



3w 



y/3{-wr/^Po{l + w) 



(3.62) 



The Universe accelerates and the scale factor blows up at finite time tc- The Universe is led to a 
big-rip singularity. 



^In the low-energy limit, for A < 1, there may exist a moment oi p = pi, for which vanishes. However, being 
different from the case oi u > —1, it does not always exist. 
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4 Anisotropy 



So far, we have investigated the expansion in terms of the scale factor a = e^. In that expansion, 
we considered also the contributions from the anisotropy c±. In this section, let us consider the 
the evolution of the shear I3±{t), and the measure of anisotropy. The shear runs as Eq. fl2.29p . 



and the measure of anisotropy is in general given by 



dVL H H [Aa3 + ^poa"^"'] 

Let us analyze this anisotropy by cases0 



(4.64) 



4.1 w > 

For w > 0, the velocity of the shear /3± has a maximum value, 

A ^Y" at ^^(^\^_ (4.65) 
A(l + w) \wKpoJ V A / ^ ^ 

At both sides of this value, /3± decays exponentially to zero. Therefore, the asymmetry of the 
spatial axes due to the shear cannot grow indefinitely at both ends. The shear does not induce any 
singular behavior in anisotropy /. However, the expansion vanishes {H = 0) at p = = X/wk 
as it was discussed in subsection 3.1. Therefore, the anisotropy / diverges at that moment. 

The scale factor approaches a constant value at this initial moment (t — oo), a « + 
j^e^ot _^ ^Yie velocity of shear becomes finite, P± ~ c-|-u;/A(l + w)a^Q = b±. Although the 

anisotropy diverges at that moment, it is not difficult to show that the spacetime is not singular 
with the metric functions a ~ aB + Ae^°^ and (}± ~ const ant. Therefore, there is no curvature 
singularity. 

For the sake of completeness, let us discuss the anisotropy at late times. At late times in the 
low-energy limit, the expansion parameter is given by Eq. (13.371) . and the measure of anisotropy 
becomes 



A 



Xa^H Xa^ I 3k 



1 {X-l)w{w + l) 
3 2A 



^2 



2 ■ 1 



1 + w 



2 



pH-^(^) } . (4.66) 



^For A < 1, there exists a bouncing moment at p = pt at which H is zero as shown in Eq. (|3.40p . At this 
moment, f3± takes a finite value because the scale factor is finite. Therefore, the measure of anisotropy / diverges. 
However, the spacetime is regular since all the metric coefficients arc finite. 
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For A > 1, the first term (cosmological constant term) dominates, and the anisotropy goes to 
zero. When the cosmological constant is absent (A = 1), as it was discussed in subsection 3.1, for 
< w < 1, the second term (linear in p) dominates and the anisotropy dies out. For w > 1, the 
third term (anisotropic term) dominates and the anisotropy becomes / — ?■ 1. 



4.2 w = 

For w = 0, the Universe approaches a de Sitter state at early times, a oc e^*, with constant 
H = a/ Acfr/3. The velocity becomes j3± ~ c±/kpq. Therefore, the anisotropy approaches a 
constant value / ~ c/npo ^Aefr/3. 

4.3 -l/3<w;<0 

For u; < 0, we observe from Eq. fl4.63p that the velocity (5± diverges as the scale factor a decreases. 
For —l<w<0, the dominant dependence is P± oc a^^ at early times. For — l/3<w<0in the 
high-energy regime (p ^ Pc), the scale factor is power-law a oc t'^/^(^+^) from Eq. (13.501) . so the 
expansion rate becomes H oc t~^. The initial anisotropy then becomes 

/ = ^oc — OC — =t— >° ^ ast-^0. 4.67 

For — 1/3 < w < 0, we observed in subsection 3.3 that there is a moment p = Pc at which H 
diverges. At that time, the spacetime becomes singular, but the anisotropy / vanishes because /3± 
is finite while H diverges. 

4.4 -l<w;<-l/3 

For — 1 < w < — 1/3 at early times, the scale factor is power-law a oc 

^-i/3w f^Qj^ (1332|1 . and 

the expansion rate is again H oc t~^. The anisotropy then becomes 
which means finite. 



4.5 W = -1 

Using the scale factor a in Eq. (I3.55P and the shear in Eq. (I3.56P , the anisotropy can be evaluated 
as 

/ = (4-69) 

cosh(V3At) 

The anisotropy goes to a constant value as t — > 0, and decays as i — >■ oo. 
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4.6 w < -1 

At early times (a the dominant dependence of the velocity of the shear is P± oc a~^. The 
scale factor is a oc (t — t^Y^^ from Eq. fl3.60p . and the expansion rate becomes if oc (t — to)^^- 
The anisotropy then becomes 

/ = ^ oc oc ) ^— 1. 4.70 

At late times (a the dominant dependence of the velocity of the shear is /3± oc a^"". The 
scale factor is a oc {tc-tf/^'^'^^'"^ from Eq. ^Mj, and the expansion rate becomes H cc (tc — t) ^ . 
The anisotropy then becomes 

(t.-r^ Hte-t)^>° ^ (ast^to). (4.71) 

5 Conclusions 

In this work, we investigated the evolution of the Universe driven by barotropic perfect fluid in 
Eddington-inspired Born-Infeld gravity. We considered both the isotropic and the anisotropic 
expansions for n > 0. 

Since EiBI gravity is the same with GR in vacuum, the evolution of the Universe at late times 
when the energy density is very low, is very similar to that in GR. For phantom matter [w < — 1), 
the energy density at late times grows, but the evolution is still similar to that in GR. 

At early times when the energy density is large, the evolution is somewhat different. For w > 0, 
the Universe starts from a "nonsingular initial state" of finite size at which the expansion rate H 
becomes zero. This was observed specifically for radiation [w = 1/3) in the original work for EiBI 
in Ref. p. 

The most interesting phenomenon arises for pressureless dust [w = 0). Even when the cos- 
mological constant is absent, A = (A — 1)/k = 0, the Universe approaches a de Sitter state at 
high-energy densities with the effective-cosmo logical constant, 

Aefi = - + (5.72) 

which provides repulsive gravity. This opens a new possibility of avoiding a singularity. At the 
final stage of collapsing dust, the high-density state may give rise to repulsive gravity, and thus 
the singularity may be not formed. This might be related with work in Ref. |2] in which the 
authors studied the interior of the compact star composed of pressureless dust. They found a 



16 



static configuration of tlie star in tlie Newtonian limit of EiBI gravity, wliicli does not exist in 
general relativity. 

The anisotropy in EiBI gravity is harmless contrary to GR. Most of the cases, the measure of 
anisotropy / dies out, or remains constant except for the initial state of the w > case. At that 
initial moment for the w > case, / is divergent, but there is no curvature singularity since the 
metric functions behave regularly. For w = 0, the initial de Sitter state has a constant value of /. 
The spacetime singularities originate mainly from the singular behavior of the scale factor a, or 
the Hubble parameter H for the w; < cases. For w > 0, the spacetime is singularity free. 
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A Energy- Momentum Conservation 

The notation for the covariant derivative used in this section is as following; 

Vj^: the covariant derivative defined by the Christoffel symbol based on g^is 
V|j: the covariant derivative defined by the connection F^^^ based on qa/s 

In order to check the energy- momentum conservation explicitly, we recall two field equations 
(fO) and ([Lll) given by 



-qq 



gf„y + kR 



fJ,U- 



(A.73) 
(A.74) 



Applying the covariant derivative on both sides of Eq. (I A. 731) leads to 



d,. 



-qq 



q q 



-gT^"" 



-q 



-gV^T^r 



(A.75) 

(A.76) 
(A.77) 
(A.78) 
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Note that in Eqs. flA.75P and flA.76|l . we used the definition of the covariant derivative for the 
tensor density |M| of weight u [TO] . 

(A.79) 



where T|^^ is the corresponding connection for the covariant derivative V^. In Eq. (]A.77P the 
following relations were used 



-g, and T^^ 



2~ ■ ■ v-9 

Next, we apply the covariant derivative on Eq. ( ]A.74|) . then we get 



--d„ 



(A.80) 



(A.81) 



where V^q^u = was used. Performing the permutation of indices in this relation (1A.81[) and 
using the definition of the covariant derivative, we get a relation. 



{1} + c, 



where C^^ is given by 

From the relation flA.82p together with flA.83p . we can express Vf^q'^" in terms of C^^p, 
Vlq^- = = d,q^^ + Tlx'' + TL'?''^ = "^iq'" + C,^V' + C,:q'^ 



(A.82) 



(A.83) 



(A.84) 
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aX 



^9 Q {CaXa ~ Cxaa ~ Caax) 



-g q t-'craA "1^9 Q ^aXa +9 q ^Xaa 



(A.85) 



Plugging Eqs. ( 1A.82I) and ( 1A.85I) into Eq. ( ]A.78p . and rearranging it by using Eq. ( ]A.81I) . we get 



^grpp.u 



(A.86) 



The quantity in the parenthesis is the Einstein tensor defined by the auxiliary metric, G'[r(g,j,)]. 
From the Bianchi identity, the left-hand side vanishes, which provides the energy-momentum 
conservation that we expected. 



(A.87) 
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B Scale Factor a{t) for p = 



We derive the scale factor a(t) explicitly for the isotropic Universe filled with dust (p = 0). Let 
us introduce a new variable, 

X2 1 X 

^ X 



z = 



z 



fB.88) 



A a/1 + Kp/A 

For the isotropic case, c± = 0, /(X, F, p) defined in Eq. f l2.25p becomes a function of z only, 

/(X, ^ /(^) = 1 + 2Az - 'iz\ (B.89) 
Therefore, Eq. f l2.24p becomes a differential equation which depends only on z^ 



dt 



where the signature ± follows from Eq. fl2.17p using Eq. f l2.23p . 

~ X ^ Y ~ 2z^ Y ~ 2z \l Qk 
Plugging Eq. (IB.SQP into Eq. (IB.QOp . we get the integral equation between z and t, 



Uz) 



dz' 



2(A - 32') 



v/2 



\z- - l)^l + 2\z' -3/^ 
The integration is performed to give 

2z 



Vl + 2A^' - 2>z' 
z'{z'^ - 



/2 



± 



{t-u). 



(B.90) 



(B.91) 



(B.92) 



Uz) 



log 



+ 



+ 



1 + A^ + VI + 2Xz - 3^2 

1 + A + (A - 3)z + V2(A- l)Vl + 2A2-3;z2 



A-1 



log- 



log- 



(l-z)[l + A + v/2(A-l)] 
A + (A + 3)^ + v/-2(A + l)Vl + 2A2-3z2 



(B.93) 



A + 1""" (1 + 2)[1-A + v/-2(A + l)] 

The arguments of the logarithm are complex in general. For A > 1, it is real if < 2; < 1. 



For < A < 1, it is real if < 2 < (A + VA^ + 3)/3. The upper limit = (A + VA2 + 3)/3 
corresponds to the minimum energy density state in which the expansion parameter becomes zero, 
H = 0. (We shall prove this later.) When this point is reached, the Universe stops expansion and 
bounces back to contract. For z > 1 and z < 0, the argument becomes complex. At z = and 
z = 1, T{z) is logarithmically divergent. For A = 1, T{z) can be written in a simpler form. 



nz) 



- 1 - tan"^ 



2z 



Vl + 2z- 3z^ 



log 



2z 



1 + 2 + ^1 + 2^-3^2 



(B.94) 
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The time evolution of z can be given by the inverse function of T as 



±J-{t-to) 



(B.95) 



although obtaining the inverse function is nontrivial. We finally get from Eq. ( IB. 911) using 
Eq. flRMj) . 

v7 



1 , 1 

- m ^ ± —= 
2 

-mz -\ — / az 
2 Qj 



fdt = - In z ± , , , , 

2 y/Q^J dz/dt 



dz 



f 



f 



z^ - 1 z{z^ - 1; ^ 

2{\-3z) l + 2\z-3z^ 
z^-1 z{z^ - 1) 



Therefore, the scale factor becomes 



a{t) 



z^l\t) 



(B.96) 



(B.97) 



[1 -^2(t)]l/3- 

Using this result, the scale factors in two limits (high- and low-energy) can be obtained with the 
aid of Eqs. ( ]B.93P and (]B.95p . which should agree with the results in subsection 3.2. 
Recasting Eq. f lB.9ip . we get 

2 



2z 



1 



— (1 + 2A^-3^2) . 



(B.98) 



When z = Zm = + VA^ + 3)/3, this equation becomes zero, i.e., tl = z/2z. By differentiating 
the scale factor a in Eq. flB.971) . one can show that H = VL = {iaXz = Zm- As z = (1 + Kp/A)"^/^ 
defined in Eq. ( IB. 881) . we have for A > 0, 

2A 



P > Pmir 



K 



(B.99) 



For < A < 1, this indicates that there is a minimum energy density of the Unverse, which is 
positive definite. (Note that for A < 1, the cosmological constant is negative.) Afterwards, the 
Universe bounces back to collapse. 

For A < 0, on the other hand, we get 

2IAI 



P 



< 



;i + a2 + |a|v/a2T3) <o. 



(B.IOO) 



Therefore, the energy density is negative definite in this case. 
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Figure 1: Plot of if vs. Q = log a for the case of w > {w = 1/3). We set c± = and po = 1- 
From the top, the lines are for A = 1.5 (dotted: positive cosmological constant), A = 1.0 (solid: 
vanishing cosmological constant), and A = 0.5 (dashed: negative cosmological constant). The 
expansion vanishes initially, H = 0. 
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Figure 2: Plot of H vs. Vt for the case of w; = 0, with the same values of parameters as in Fig. 1. 
H becomes constant at early times, which indicates that the Universe is in the de Sitter state. 




Figure 3: Plot of H vs. Vt for the case of —1/3 < w < (w = —1/6), with the same values of 
parameters as in Fig. 1. There are two universes split by the singular point p = pc aX which H 
diverges. 
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Figure 4: Plot of H vs. Q for the case of — 1 < < —1/3 {w = —0.6), with the same values of 
parameters as in Fig. 1. 




-1.0 -0.5 0.5 1.0 1.5 2.0 

Figure 5: Plot of H vs. fl for the case oi w < —1 {w = —1.5), with the same values of parameters 
as in Fig. 1. 



24 



